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1. Introduction 

Despite the fact that only very few successful attempts to produce hole-doped crystalline 
fuUerene have been reported (see Panich et al. [T1I21IS]), it is interesting to investigate the 
properties of such compounds. Their physical behavior will to large parts be determined 
by the ground state properties and the low energy excitations of positively charged 
molecular ions. Besides that, fullerene molecular ions can be realized in solution, and 
to some extent also in vacuum. 

The physics of charged fullerenes is determined by their high orbital degeneracy 
(3 in case of electron doping and 5 in case of hole doping for Ceo), and the resulting 
interplay between the Jahn- Teller (JT) effect and Hund's rule. The latter favors high- 
spin ground states while the former splits the orbital degeneracy by a JT distortion of 
the molecule, leading to a low-spin ground state. This in turn will influence the tendency 
of hole-doped systems to become e.g. insulating or magnetic or superconducting. 

This interplay can be observed in the electron-doped compounds, some of which 
show metallic behavior and can become superconducting, such as A3C60, (A=K, Rb, 
Cs) jHE], others are insulating like, e.g. Na2C6o |Hj or A4C60 0- Recent data on the 
charge transfer compounds (AsF6)2C6o and (SbF6)C6o, which contain nominal Cgo^'*' 
seem insulating and compatible with a magnetic ground state of the Ceo^"*" ions jSj. 

In this paper we report a theoretical study of the multiplet structures of isolated 
Cgo"^ ions. Due to the open shell nature of the charged molecules, the multiplet energies 
cannot be obtained reliably by more standard methods such as ROHF [Restricted Open- 
Shell Hartree-Fock) or LDA [Local Density Approximation ) calculations. We apply ab 
initio CAS SCF [Complete Active Space Self- Consistent Field) calculations, which are 
the appropriate generalization of the Hartree-Fock (HF) method to orbitally degenerate 
systems jSJ. 

We compare our results with previous studies, which were based on the semi- 
empirical calculations OE], model-SCF calculations ^1], limited CI [Configuration 
Interaction) calculations fl2\ IT3] and on constrained LDA calculations p3] . 

2. Model Hamiltonian 

As the starting point for describing the low energy physics of doped fullerenes one needs 
a model Hamiltonian which describes the interactions among the tiu electrons in case 
of electron-doping and or the hu holes in case of hole-doping, coupled to the vibrations 
of atoms in the molecule. Such a model was recently proposed in Ref. ^3] and shall 
only be briefly introduced here. 

The model Hamiltonian for a single molecule reads 

H = Hq + ifvib + -^e-vib + -^e-e, (l) 

where 

Ho =eY,cl^c^m, (2) 
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He-vih = ^ ^ 2 Cj^jI^^, QiA/iClm'^am' , (4) 
amm' riAfi 

He-e =^Y.Y. ^<ry{,m, m'; n, n') c),„cj^,„, c^'n'Can, (5) 



cr,cr mm' 



are respectively the single-particle Hamiltonian, the molecular vibration energy in the 
harmonic approximation, the electron-vibration coupling in the linear JT approximation 
[inilini, and the mutual Coulomb repulsion between the electrons. The denote the 
creation operators of either a hole in the HOMO [Highest Occupied Molecular Orbital) 
or an electron in the tiu LUMO [Lowest Occupied Molecular Orbital), described by the 
single-particle wave function (/^^^-(r). The index a is the spin projection, m and n label 
the component within the degenerate electronic HOMO/LUMO multiplet, and i counts 
the vibration modes of symmetry A (2 Ag, 6 Gg and 8 modes). C^'^ are Clebsch- 
Gordan coefficients of the icosahedral group, for coupling hu (holes) or ti„ (electrons) 
states to vibrations of symmetry A. The index r =1,2 is a multiplicity label, relevant 
for Hg modes only p!5t ITT] . QiAfi and PiA/i are the molecular vibration coordinates and 
conjugate momenta. 

The Coulomb matrix elements are defined by: 

Wa,a'[m,m -,71,71') = Jd^r Jd^r' (p*^^[r) (p*^,^,[r') Ua,a'[r,r') ipna[r) (Pn'a'[r'), (6) 

where Ucr,a'[r,r') is a generally screened Coulomb interaction. In Ref. |21 it was 
demonstrated that this interaction can, without loss of generality in the given subspace, 
be expressed in terms of a minimal number of physical parameters as: 

rr'A \ M / 

where we defined the operators: 



w 



rAfi ^ f^rA/i A (o) 



The single-particle term A, which results from rearranging the field operators, is a 
function of the Coulomb parameters defined below. 

For electrons, the allowed Coulomb parameters, which generalize Slater's Coulomb 
parameters for atoms |18j, arise from the product ti^^^iM = ^g®Hg and we define thus 
Fi = and F2 = F^o . For holes, the relevant product is hu®hu = Ag®Gg® 2Hg, 
giving rise to five parameters Fi = F^3,F2 = F'^^F^ = F^'^'^o,F4 = F^'^-^s and the 
cross-term F5 = F^'^'^s. 

Interaction parameters with a more transparent meaning can be further introduced 
by looking at the averaged energies of a multiplet of given charge n and spin S. We find 
that 

Tiin — 1 ) 

E^,,[n) =Eo + en + U \ ^ (9) 
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Table 1. Multiplet energies of the N electron states (tiu)^ relative to the average 
energy eN + UN{N + l)/2. 



N 


S(S+l)^ 


£'inult 


N 




£'inult 


2,4 


'A, 


4J 


3 




2J 






J 













-J 






-3J 



and 

E,.e(ri, 5) - E,,,{n, S - 1) = 2JS (10) 

for both electrons and holes. To make contact with standard notation established 
historically, we rescale the exchange parameter for electrons as J = 4 J/5, but keep 
J = J for holes. It should be noted that here the on-site Coulomb repulsion parameter 
U is defined with respect to the multiplet averaged energies, while traditionally the 
"Hubbard" U (we call it Uq in this work) is defined with respect to the ground-state 
energies of each charged system, i.e. 

Uo = I{n) - A{n) = I{n) - I{n + l) = EQ{n + l) + E^{n-l) - 2Eo{n), (11) 

where I{n) and A{n) are respectively the ionization potential and the electron affinity 
of the n-charge ion. 

The parameters introduced above can be expressed in terms of the original 
Coulomb-parameters, as: 

Utiu = ^ (^1 - ^2) (12) 

■Jtiu = 2^2 (13) 

1. 
-( 
2 



= lu + 2J (14) 



for electrons, and 

u, =!i_^Ii^!a_E± (15) 
" 5 45 9 9 ^ ' 

Jku =\f2 + ^^{F, + F,) (16) 

^^^^ = \u + \j (17) 

for holes. Analytic expressions for the eigenenergies of Hamiltonian Q for holes were 
partly already reported in Ref. ^3]. For completeness, we provide the eigenenergies 
for the N electron states in Table d Our scope in the rest of this paper will be the 
calculation of all these Coulomb parameters. 

3. Computational details 

In order to find the multiplet splittings we performed CAS SCF |8j calculations using 
the quantum chemistry package GAMESS [19]. If not stated otherwise, all calculations 
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Figure 1. The orbital degeneracies and splittings close to the HOMO-LUMO gap. 



LUM0+l(Ti3) 
LUMO (Ti„) 



1.69 eV 



7.54 eV 



HOMO (H„) 
HOMO-1 (Eg) 



1.67 eV 



presented here have been performed using the 6-31 Gaussian basis set [6s3splsp]/[ls2sp] 



The energetic picture of the canonical orbitals, shown in Figure ^ motivates the 
choice of the active space. For the positive ions, we consider only states with holes 
delocalized in the 5-fold degenerate HOMO shell, and for the negative ions, the states 
with electrons in the 3-fold degenerate LUMO shell. This is justified by the large HOMO- 
LUMO gap of the isolated Ceo molecule 7.54 eV. The separation between the HOMO 
shell and the HOMO-1 shell is 1.67 eV. A similar gap of 1.69 eV is found between the 
LUMO and LUMO-1 shells. We checked moreover that including the HOMO-1 and 
LUMO-l-1 shells respectively did not change the results for one and two holes in the 
HOMO and one and two electrons in the LUMO. Thus, all results are obtained with 
only one degenerate shell in the active space. The natural orbitals were optimized with 
respect to the averaged energy of all possible multiplet states regardless of the spin 
multiplicity. The geometry for the neutral Ceo was optimized with the RHF {Restricted 
Hartree-Fock) method with perfect icosahedral symmetry. The hexagon and pentagon 
bond lengths obtained are 1.474 A and 1.391 A respectively. These are not far from 
the experimental crystalline values obtained from the neutron powder diffraction study 



(1.455 A and 1.391 A f^P), and NMR measurement of C-C bonds (1.450 A and 

1.40 A jS!). II 



II In order to further improve these results, one would need a bigger basis-set and inclusion of dynamical 
correlations in the calculations. The importance of these effects on the geometry and vibrational modes 
of Cgo is discussed in Ref. \2'6\ . Performing dynamically correlated calculations within the perturbation 
theory, for instance on the MP2 level {second order Mdller-Plesset) |24[ I25| . would be computationally 
very expensive for the systems studied here. MP2 is a Rayleigh-Schrodinger perturbative method which 
uses the Hartree-Fock nonperturbed wave function and Hamiltonian as a zeroth order approximation. 
The perturbed functions depend on single-particle energies therefore it is a so-called "method with 
dynamical correlations" . The semiempirical-approaches by the use of spectroscopic data for the 
estimation of two-electron integrals and the LDA method by the parametrization on the QMC data 
{Quantum Monte Carlo) also partially include the dynamical correlations. To estimate the effect of 
dynamical correlations, we report some results where we include them partially for the Cqq ion. 
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Although the JT effect is crucial for the physics of Cgo ions ^Hj, these changes 
of the geometry are relatively small and their influence on the Coulomb integrals (see 
Eq. (jni)) should be fairly small. In the present calculations we ignore therefore those 
distortions and fix the geometry as frozen to that of neutral Cgo- 

4. Results 

First we report the ground state energies for the Cgo ions relative to the ground state of 
the neutral system. They are shown in Table El together with the ionization potentials 
and, if known, the corresponding experimental values. We report also the Hubbard t/'s, 
calculated from the ground state energies. 

The calculated ionization potential (IP) in the +1 charged case 10.93 eV is in a 
good agreement with the experimental value 11.46±0.05 eV jSHj. Also for the neutral 
molecule the calculated IP at 7.88 eV is close to the measured 7.6±0.2 eV [21] . Not 
surprisingly, the IP for the negative ion at 0.63 eV differs much from the experimental 
2.65±0.05 eV The negatively charged ion is theoretically much more demanding 
since the additional electron moves in a less binding potential, is very delocalized and 
very sensitive to correlations with the other electrons. Therefore the basis set and the 
relaxation effects are more important in negative ions. To check the effects of the basis 
set and of additional correlations, we also performed calculations with the minimal basis 
set of order 3-21 ([6s3p]/[2slp]) jSH], which also allowed for a R0MP2 {Restricted Open 
Shell MP2) calculation, with frozen 175 occupied orbitals and 10+1 correlated electrons. 

Table 2. Ground state energies (in eV), ionization potentials and Hubbard J7's (f/o, 
obtained from the ground state total energies according to eq. for the charged 

fuUerene. 



Charge 


^GS 


/ 






-6 


43.157 


-14.974 






-5 


28.183 


-12.009 




2.97 


-4 


16.174 


-8.999 




3.01 


-3 


7.175 


-5.415 




3.58 


-2 


1.760 


-2.393 




3.02 


-1 


-0.633 


0.633 


2.65 ± 0.05 


3.03 





0.000 


7.883 


7.6 ± 0.2 




1 


7.883 


10.925 


11.46 ± 0.05 


3.04 


2 


18.808 


13.989 




3.06 


3 


32.797 


17.052 




3.06 


4 


49.849 


20.110 




3.06 


5 


69.959 


24.336 




4.26 


6 


94.295 


27.403 




3.07 


7 


121.698 


30.492 




3.09 


8 


152.190 


33.591 




3.10 


9 


185.781 


36.6728 




3.08 


10 


222.453 
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The restricted basis set reduced the IP of Cqq, as expected, to -0.03 eV. Inclusion of the 
dynamical correlation increased the value from -0.03 eV to +0.43 eV. We checked that 
inclusion of 5- fold degenerate orbitals on HOMO-1 level in the active space does not 
change the IP for the CgQ ion by more than 0.06 eV. Also the inclusion of higher 3-fold 
degenerate orbitals gives minor changes. The excitations from these orbitals contribute 
to the multiplets of negative ion with total weight smaller than 2% and do not change any 
of the splittings. The better agreement with experiment obtained by other authors by 
means of semiempirical approaches [01 EDI EI] for CgQ could be partly due to the fact that 
the parameters used there are derived from spectroscopic data. CgQ ions were observed 
in gas phase measurements [S21 EB] and found to be stable in semiempirical calculations 
jniEni- However the experimental results were questioned by other group [33]. In print 
of principle, our calculations are not ideally suited for negative ions. First of all, the 
weakly bound electron state is delicately sensitive to polarization and to correlations, not 
very accurate in DPT. Secondly, the basis set becomes critical unless it is extended very 
much, which is beyond our scopes. Earlier calculations by Razafinjanahary et al. jHS] 
are probably of better quality in that respect. A comparison to the previous calculations 
by means of DPT {Density Functional Theory) given in Table 4 of Ref. [SS] supports the 
importance of electronic correlations. The IPs collected in that paper strongly depend 
on the exchange-correlation functional as one can see from values: 2.0 eV jHEl; 1-9 eV 
[SS], 2.7 eV inni and 2.8 eV [SZI for CgQ and values: -1.3 eV [36j, -1.2 eV [SS], -0.4 eV 
and -0.3 eV [BZl for Cgg obtained respectively with functional: Xa, VWN (Vosko- 
Wilk-Nusair) , BH (Barth-Hedin) and PZ (Perdew-Zunger) . Our IP of CgQ found to be 
-5.42 eV should be compared with value -4.4 eV obtained with the VWN functional 
from the modified MT [Muffin-Tin Approximation) calculations by Razafinjanahary et 
al. [S3], t 

The bare molecular Hubbard f/'s (Uq) are all of the order of 3.01-3.10 eV for both 
electrons and holes (see Table 2). Exceptions are the 3-fold negative and the 5-fold 
positive systems. Here the ?7's are considerably larger. The parameters Uq obtained for 
the charge n were derived from the ground state energies of molecules with charge n — 1 
and n + 1. Due to the particle-hole symmetry, the multiplet splittings of ions with charge 
one more and one less than the charge at half-filling are identical (the energy gaps on 
right and left from half- filling are equal). This is not the case of any other charge, where 
the "neighboring" multiplet splittings are different and the ground energies of the charge 
multiplets which are more rich in the number of states, are lowered in comparison to 
the ground energies of multiplets with smaller number of states. The aforementioned 
asymmetry leads to smaller values of Uq away from exact half filling. 

^ The question was raised in Ref. |35j about possible electron states trapped inside the Cgo molecule. 
If such states existed, a calculation like the present where all basis functions are carbon-centered, and 
none is molecule-centered, could very well miss such a trapped state. However, the trapped state would 
not be missed in a plane-wave (PW) calculation. Direct comparison of our tiu LUMO and t2g LUMO-l-1 
derived states with those of the PW calculations of Cgo in Ref. [2HII31I shows perfect agreement, with 
no additional low laying states in the PW case that could be taken as trapped states (states with the 
spherical symmetry). 
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Table 3. Multiplet energies (in Merv) for 2,3 and 4 electrons in the LUMPY, 
and energies reconstructed from the model (eq. (11)-(13)) using the charge-averaged 
parameters eq. (9,10). 



A 



N=2 



N=4 Model 



A 



N=3 Model 



^Tig -114.1 



^Hg 114.1 
^Ag 456.6 450.9 



112.7 -113.5 
112.7 113.5 
453.9 



U 



-340.7 -340.4 
0.0 0.0 
227.1 226.9 



Table 4. Multiplet energies (in meV) for 2,3,7 and 8 holes in the (/i„) HOMO, 
and energies reconstructed from the model (eq. (14)-(16)) using the charge-averaged 
parameters eq. (9,10). 



A 


N=2 


N=8 


Model 


A 


N=3 


N=7 


Model 




-134.8 


-125.2 


-133.8 


J-bu 


-388.9 


-382.3 


-390.3 


J^ag 


-127.9 


-139.7 


-128.7 


^ au 


-385.8 


-394.8 


-385.2 


^Gg 


-122.4 


-122.5 


-123.8 




-378.8 


-382.0 


-380.2 


Gg 


52.0 


63.6 


57.7 


au 


-113.6 


-99.2 


-109.1 


'Hg 


57.4 


67.1 


62.2 




-110.5 


-107.8 


-104.0 


'Hg 


463.3 


448.4 


455.9 




-105.6 


-95.3 


-101.3 


'Ag 


1022.3 


1029.7 


1026.0 




-102.2 


-91.3 


-97.4 










^Gu 


191.2 


189.1 


190.1 










J-lu, J2u 


199.7 


196.0 


197.7 










Hu 


398.1 


386.9 


392.1 










'^Gu 


399.1 


388.1 


393.1 












767.9 


772.1 


768.6 



The key results of the above described CAS SCF calculations are the multiplet 
energies of the fullerene ions. They are reported in Table 01 for electrons and Tables |3] 
and El for holes. Since GAMESS does not support symmetries as high as icosahedral, 
the assignment of the symmetry labels to the states in CAS SCF calculations was done 
indirectly, using the degeneracies of the states. In this way, the three-fold degenerate Ti 
and T2 representations cannot be distinguished. However, due to the analytical results 
of the model Hamiltonian, reported in Ref. pH], they can be assigned once a certain 
sign of the parameter F5 is chosen. Hence, in these cases we label the states by Ta/h 
where a/h represents 1 or 2 respectively. In the negative ions, no such ambiguity arises. 

We find that the spread of the multiplets is about 0.6 eV for all negatively charged 
states, 1.2 eV for 2, 3, 7 and 8 holes, and 2.3 eV for 4-6 holes. Comparing the spectra for 
N and 6-A^ electrons, or N and 10- holes respectively, one finds that the particle-hole 
symmetry in the occupation of the degenerate orbitals is approximately preserved in 
both cases. 

The semi-empirical calculations for multiplets of Cg^J" ions and for Cgo^ and Cqq ions 

Clearly however, this approximate symmetry is a result of out single-orbital approximation, and would 
disappear once all other levels would be included. In this sense, it is mostly of academic significance. 
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Table 5. Multiplet energies (in meV) for 2,3,7 and 8 holes in the HOMO, 
and energies reconstructed from the model (eq. (14)-(16)) using the charge-averaged 
parameters eq. (9,10). 



A 


N=4 


N=6 


Model 


A 


N=5 


Model 




-768.3 


-770.1 


-769.5 




-1282.7 


-1282.6 




-395.7 


-389.8 


-393.1 


Hu 


-712.6 


-712.4 


3 TT 

Hg 


-394.7 


-388.9 


-392.2 




-331.3 


-331.1 




-209.5 


-199.7 


-204.9 




-330.8 


-330.6 


1 A 


-205.5 


-196.2 


-201.2 




-318.7 


-318.7 


Lfg 


-197.0 


-195.2 


-196.2 




-314.2 


-314.2 


Tbg 


-194.9 


-196.3 


-195.6 


2 A 


-139.1 


-139.0 


1 


-193.9 


-189.7 


-192.1 




-132.8 


-132.8 


1 TT 

Hg 


86.2 


88.9 


87.7 




-131.5 


-131.4 


-Lag 


92.0 


91.0 


91.7 


4^^ 
a u 


-127.7 


-127.8 


^ TT 

Hg 


93.8 


97.9 


95.5 


4t~' 
^ lu 


-122.7 


-122.7 


\rr 

J-hg 


93.5 


97.2 


96.7 


2 TT 


143.6 


143.8 


1 


103.5 


100.6 


100.6 


2 TT 


156.1 


156.1 


S TT 

Hg 


103.3 


98.3 


101.9 


J-hu 


158.9 


158.9 


Tbg 


104.0 


100.6 


102.4 


2 TT 


159.0 


159.0 


3 TT 

Hg 


105.9 


100.6 


103.4 


^ au 


162.9 


162.8 


n 

y 


382.6 


382.0 


383.4 


^Tau 


164.3 


164.4 




384.3 


383.8 


383.4 


^Hu 


166.9 


166.9 


^ag 


385.1 


384.7 


385.0 


^Tbu 


167.9 


167.9 


'Tbg 


385.6 


387.9 


386.7 


'Tbu 


443.9 


443.7 


'Gg 


392.3 


388.8 


390.9 


J- au 


453.4 


453.4 


[Hg 


564.3 


569.9 


567.2 




455.7 


455.6 


Gg 


578.4 


581.1 


580.2 


^Gu 


460.1 


459.9 


[Hg 


590.4 


583.1 


587.2 




463.9 


463.6 


Gg 


591.3 


585.2 


588.3 


^Gu 


646.8 


646.7 


'Hg 


970.3 


966.9 


969.0 




651.3 


651.1 


'Ag 


1536.7 


1540.7 


1539.2 




1026.3 


1026.2 



give larger splittings than ours probably due to unequal occupation of degenerate 
orbitals. The limited CI calculations with HF reference orbitals do not take care of the 
right degeneracies due to the splitting of ti^ shell when it is partially occupied by the 
electrons. The discussion of this fact is given in Ref. [12]. The CAS SCF calculations 
are free of symmetry related problems. Also, excessively large splittings were obtained 
with the restricted configuration interaction constructed using the INDO {Intermediate 
Neglect of the Diatomic Overlap, for the method see jUIj) orbitals [Tllj . 

Recently a complete analysis of the multiplets of all Ceo ions was presented [TT], 
where the multiplet energies were evaluated from a multipole expansion of the Coulomb 
interaction. The Coulomb integrals were calculated with respect to symmetry-adapted 
model wave functions. The obtained multiplet spectra are in good agreement with our 
results, apart from a larger T1/T2 splitting, and a change in the level ordering. We will 
comment on this point later, when discussing the Coulomb parameters extracted from 
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Figure 2. Average multiplet energies and the fitted parabola for the (tiu)^ states 
for electrons and the [hy)^ states for holes. 
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the present results. 

We observe a degeneracy between the ^Ti„ and the ^T2„ states of the three hole 
multiplets. This apparently accidental degeneracy is indeed well studied in the literature 
im 1321 Hn] and attributed to the permutational symmetry. 

5. Determination of the model parameters 

According to Eq. 0, the multiplet-averaged energies of the model are described by the 
parabola: 

U 71^ 

E^,,{n) =Eo + {e- -)n + U—, (18) 

where n denotes the number of electrons or holes respectively. 

Figure El displays the calculated average energies, together with the fitted parabolas: 

El~^{n) = - 61393.320 - 2.145n + 1.557n^ (19) 
E^^^{n) = - 61393.268 + 6.245n + 1.599^2. (20) 

From these parabolas the average Coulomb repulsion ?7ave and the single particle energies 
e are obtained. These are displayed in Table IHl 

This Coulomb repulsion f/ave, of about 3 eV for both electrons and holes, is in good 
agreement with previous calculations |^ IT^ . The single-particle energies, extracted 
from the parabolas, are very close to the single-particle levels of the neutral system. It 
should be noted that this is not an obvious finding, since the single-particle energies, 
extracted from the fit, derive from a calculation of the charged molecules. In the model 
Hamiltonian, the first ionization potential is equal exactly to the single particle energy 
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Table 6. Coulomb parameters end Koopmans single particle energies for electrons 
and holes (in eV). 





UavG 


e 


crhf 


electrons 


3.1141 


-0.5881 


-0.4898 


holes 


3.1928 


-7.8477 


-8.0270 



Table 7. Exchange splitting J of the electron states (iiu)^, for each charge N, as well 
as it's averaged value (ave). The values labeled by limited CI are extracted from 
INDO-I-CI from [TUI, "direct" refers to Nikolaev et.al. mi, LDA to Liiders et al. [H] . 
MNDO-I-CI to (301 ■ AH parameters are given in meV. 



Method 


N=2 


3 


4 


ave 


CAS SCF 


114.1 


113.6 


112.7 


113.5 


limited CI 


29-146 


8-388 


8.5-380 


15-305 


INDO-fCI 


120 


185 




153 


"direct" 








99, 95 


LDA 








32 


MNDO-I-CI 


50 






50 



of the holes. In fact, Koopmans theorem states that the HOMO level gives a good 
approximation for the ionization potential. In general Koopmans theorem does not 
perform very well for the electron affinity, due to a bad description of the virtual states 
by means of the HF method (see ^2] P- 127-128). Calculations involving charged states, 
in which the LUMO shell of the neutral molecule is actually occupied, should give much 
better values. Here we find, however, that the single particle energy of the electrons, 
determined from the fit, which represents the electron affinity of the model, is very close 
to the LUMO level obtained by means of the RHF method. 

Next, for the N electron states {tiu)^ only the exchange parameter J remains 
to be determined. Knowing the analytical eigenenergies (Table Q), it can be directly 
extracted from the calculated multiplets given in Table El The exchange constants Js, 
obtained in this way, are given in Table [71 We compare them to other results by means 
of the LDA method, and to model calculations ("direct" in Table j?!) and to the limited 
CI method, and finally to the small configuration interaction performed on the INDO 
orbitals. Our results are close to those obtained by the model calculations |TT] , although 
slightly larger. By contrast the values of Hund's rule exchange J obtained by LDA is a 
factor of 3 smaller The reason for the large discrepancy with the LDA calculations 
is not clear. * 

The J parameters shown in Table 13 for limited CI are extracted according to 
Tabled from the multiplet splittings of either the lowest states or of the whole multiplet 

* This could be due to correlation effects included in the effective Kohn-Sham potential or due to the 
symmetry or self-interaction problems present in LDA calculations. 
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within the ti^ shell. These CI parameters depend strongly on the charge and the 
choice of states used for the derivation. This is due to the fact that the configuration 
interaction approach does not treat correctly the symmetry of ti„ shell, removing the 
3-fold degeneracy when the shell is partially occupied by the electrons j^. CAS SCF 
calculations give correct degeneracies and the obtained J parameters are independent 
on the charge and the choice of states (due to the multiplet structure showed in Table HJ. 
The fact that limited CI calculations underestimate splittings of the lowest states and 
overestimate splitting of higher states is not only due to the broken symmetry, but also 
due to the use of HF single-particle orbitals optimized for the ground state. Contrary 
to that, in our CAS SCF calculations we optimized the molecular orbitals in such a way 
that the total energy averaged over all states within the multiplet was minimal. 

Other authors, by means of the MNDO+CI method {Modified Neglect of Diatomic 
Overlap, for the method see ^Sj), derived J = 30 meV for the neutral Cgo molecule and 
J=50 meV from the Cqq triplet-singlet splitting j^HI- The ab initio SCF calculations 
point to a J value of order 100 meV [IHj. The parameters obtained from the MNDO-fCI 
method are smaller than the limited CI results probably due to the use of semi-empirical 
orbitals (direct fit on spectroscopic data) for the CI construction. 

For the positively charged states, the situation is both more complicated, and less 
explored. Due to the uncertainty of the Ti, T2 assignment and the multiple occurrence 
of states of the same symmetry but different energies, the parameters cannot be fitted 
directly to the multiplet energies. We therefore refer to the average energies for a given 
symmetry, shown in Table |H1 As can be seen from this table, the averaging over the Ti 
and the T2 states depends only on F2, F3 and F4, while F5 determines Ti - T2 splitting 
only. The parameters F2 - F4 can be determined completely, while F5 can only be found 
up to a sign. As it turns out, F5 is however small enough to make this uncertainty 
irrelevant. 

We fit the multiplet energies for each charge separately, allowing for charge- 
dependent model parameters, which are given in Table IHl The parameters F2,F^ and 
F4 show a weak monotonic dependence on charge. This is in contrast to the result of 
Ref. JT] whose parameters are basically determined at A^=0, and then used for all N. 
From this perspective it becomes clearer why the Tiu/T2u splittings of Ref. are 
much larger than ours. 

The substantial difference of our CAS SCF results with the LDA parameters of 
Ref. ^M, for hole-doped ions, also shown in Table El is similar to the discrepancy noted 
in the electron-doped case. 

6. Summary 

The central aim of this work is a calculation of all Coulomb parameters and multiplet 
energies for C^g" io'^s including a proper treatment of the orbital degeneracies, within 
the icosahedral symmetry of Ceo- This has been done within the CAS SCF framework. 
The model parameters for positively charged Cgo ions are substantially larger than 
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Table 8. The Coulomb multiplets for C^o+, averaged over states of the same 
representation, as a function of the e-e parameters. The model Hamiltonian obeys 
particle-hole symmetry: therefore the multiplet energies for > 5 holes equal those 
for (10 — iV) holes. The non particle-hole symmetric contribution [eN + UN{N — l)/2] 
is left out in this table. 
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Table 9. The model parameters for the positively charged states, determined from 
the multiplet spectra. We also give a comparison to the LDA calculations of Ref. ^1] 
and the model (direct) calculations of Ref. All parameters are given in meV. 





N=2 


3 


4 


CAS SCF 
5 6 


7 


8 


ave 


direct 
II III 


LDA 




188.8 


190.3 


191.9 


193.5 


195.1 


196.6 


198.2 


193.5 


193 


185 


105 


Fs 


591.1 


589.0 


586.1 


584.3 


581.3 


580.2 


577.1 


584.2 


640 


625 


155 


Fi 


178.0 


180.1 


182.3 


184.4 


186.8 


188.7 


191.1 


184.5 


180 


172 


47 


\F,\ 


3.1 


1.4 


0.7 


2.2 


4.3 


5.6 


7.6 


2.3* 


13 


13 





J 


191.7 


191.9 


192.1 


192.4 


192.5 


193.0 


193.1 


192.4 


203 


197 


60 



the LDA results as in the negatively charged molecules. Our results on the other 
hand, are rather closer to model studies by Nikolaev jTTj. A possible reason is that 
both approaches do not take into account the dynamical correlations, which are very 
important, as was shown in Ref. [22j. In the LDA and MNDO methods, where some 
dynamical correlations are included, the multiplet splittings are much reduced. The 
limited CI method led in previous calculations jT^] to much smaller or much larger J 
values depending on the choice of states and on the charge case used for the derivation of 
parameters. This could be due to the artificially broken symmetries when the electrons 
occupy partially the molecular orbitals. The parameters obtained from the CAS SCF 
method, involving a correct treatment of symmetry, do not depend on the charge or 
choice of states within the multiplet used for the derivation. 

The mean Coulomb parameters, U, obtained in this work from averaged multiplets, 
are for positively charged fullerenes about 3.2 eV. Similarly for the negative ions the 
obtained U parameters are about 3.1 eV. The charge specific Hubbard Uq are generally 
close to the average, except at half filling where Uq is not influenced by J. Here ?7o=3.58 
for C^o and Uo=4.26 for C^^. 

The exchange parameters, J, for the negative ions are somewhat basis set dependent 
due to the localization effects discussed earlier, and their values are of about 114 meV 
for 6-31 basis set and 146 meV for the minimal basis set |22j. The more extended 
basis set contains more expanded Gaussians, while the smaller basis set constraining 
the electrons to a smaller volume, is likely to increase the exchange parameter J. The 
exchange parameters for hole-doped molecules are less basis set dependent: they are of 
order 192 meV for the 6-31 basis set and 208 meV for the minimal basis set. These 
values probably still represent an approximation in excess and it seems possible that a 
further slight decrease could be found with substantially larger basis sets. ^ 

Experimentally, the multiplet structure of Ceo ions in vacuum is still, as far as 
we know, inaccessible. In particular probing the magnetic state of ions in a beam 

tt The quality of Gaussian basis set type is more dependent on the well chosen exponents than on their 
number. There is no possibility for the saturation of a basis set similar to that can obtained with the 
plane- wave methods. 
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appears to be difficult. Probing charged fullerene ions embedded in a rare gas matrix 
would seem a better possibility, that could be explored with the purpose identifying 
the ground states properties for the various charged species, possibly with EPR ff. The 
Coulomb couplings of Cqq ions play an important role in the solid state compounds. In 
particular in alkali fullerides the Cgo molecule is negatively ionized while in Lewis acid 
fullerene salts ^1 12] it should be positively ionized. NMR longitudinal relaxation time 
data of basically all alkali fullerides indicate a " spin gap" , separating the low spin Cqo 
ion ground state from high spin excitations jUj. These spin gaps are typically 100-150 
meV, compatible with a value of J of e.g. 50-70 meV. Comparison with our calculated 
value of 113 meV indicates an overestimate. While that may be due to a genuine solid 
state effect, we cannot rule out the possibility that the agreement could be improved by 
further extending the used basis set. 

Finally, our predictions of a large J for positive fullerene ions, which also imply 
that they should be with great probabihty magnetic [H] should be directly accessible 
by further experiments on the acceptor salts Pfl Oj as well as in the recently discovered 
C^Q and Cgo^ in solution |1H] . 
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